Excitation of the Dissipationless Higgs Mode in a Fermionic Condensate 
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The amplitude mode of a fermionic superfluid, analogous to the Higgs Boson, becomes undamped 
in the strong coupling regime when its frequency is pushed inside the BCS energy gap. We argue that 
this is the case in cold gases due to the energy dispersion and nonlocality of the pairing interaction, 
and propose to use the Feshbach resonance regime for parametric excitation of this mode. The results 
presented for the BCS pairing dynamics indicate that even weak dispersion suppresses dephasing 
and gives rise to persistent oscillations. The frequency of oscillations extracted from our simulation 
of the BCS dynamics agrees with the prediction of the many-body theory. 



The observation of resonance superfluidity in cold 
atomic Fermi gases P, 0] at magnetically tunable Fes- 
hbach resonances Q opened new avenue of exploring the 
many-body phenomena. Similar to the earlier work on 
cold Bose gases which triggered studies of fascinating col- 
lective phenomen a [a, fermionic pairing at Feshbach 
resonances 0, H S El El E2] presents new opportunities. 
In particular, the high degree of coherence of trapped 
atoms, and the possibility to control particle interaction 
in situ on the shortest collective time-scale, the inverse 
Fermi energy [J], can facilitate exploring new regimes 
which arc difficult to realize in solid state systems. 

The theory of fermionic pairing predicts two princi- 
pal collective modes intrinsic to the condensed state. 
One is the massless Bogoliubov-Anderson mode related 
to the order parameter phase dynamics. Being a Gold- 
stone mode, it manifests itself in hydrodynamics in the 
same way as in Bose systems, and was probed recently 
in the experiments on gas expansion and oscillation in 
In addition, there exists a second 
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fundamental elementary excitation 17|, |l8|, [l9| , related to 
the dynamics of the order parameter modulus |A|. No- 
tably, this excitation is unique to fermionic pairing and 
has no counterpart in Bose systems }19j . This massive 
excitation, characterized by a finite frequency, is anal- 
ogous to the Higgs Boson in particle physics. Like the 
latter it remained elusive, for a long time evading direct 
probes, although some indirect manifestations have been 



discussed [2JJ, l21| . The main obstacle to the detection 
of the Higgs mode in superconductors is that it is es- 
sentially decoupled from the phase mode responsible for 
hydrodynamics and superfluidity 

In this work we propose to use the dynamical con- 
trol of pairing interaction demonstrated in Refs.[l|, 0] for 
parametric excitation of the Higgs mode. We argue that 
fermion superfluidity in the strong coupling regime real- 
ized near Feshbach resonance represents a distinct ad- 
vantage, since in this case the Higgs mode is pushed 
inside the superconducting gap, htu < 2A, which elim- 
inates damping due to coupling to quasiparticles. We 
demonstrate that this mode can be excited by a time- 
dependent pairing interaction, as illustrated in Fig.[TJ In 



0.05 



a) 
5 -0.05 



*0.05 



g; o 

CD 
C 
LU 



0.05 



11 


III 


(b) 






w w ▼ w w w 



w 

c 

i 

n' 



CD 
T3 
13 

E 

CO 
O) 

£Z 

cc 
Q. 



0.5 





irA ft ft- 


A(c)A A A A A A 


A + 1\ l\ 




\l a\J v 


Time (A" F 1 ) 





w 



•0.5 







10 



20 



30 



40 50 



FIG. 1: Non-decaying Higgs mode in a Fermi gas with 
energy-dependent pairing interaction excited by the interac- 
tion switching from gi at t < to g at t > 0. Shown are 
the time and energy dependence of the pairing amplitude (a) , 
the a;-component of the pseudospin vector (b) , and the pairing 
amplitude at the Fermi energy (c) as obtained from the model 
(El),© at g = 0.43, g t = 0.23, an = a 2 = 0.5, j/W = 0.01, 
Af/W — 0.016. Note the initial transient of few periods, ex- 
hibiting some dephasing in pseudospin dynamics (b) , followed 
by synchronized collective oscillations of fermion states. 



contrast, the BCS theory at weak coupling predicts the 
Higgs mode frequency right at the edge of the quasiparti- 
cle continuum, Tiuj = 2 A 1911 . which leads to collisionless 
damping of this mode [22|, |28( . 

The departure from the behavior at weak coupling 
arises from the change in the character of pairing inter- 
action in the strong coupling regime, in particular due 
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to its finite spatial radius and frequency dispersion. Spa- 
tial nonlocality of pairing interaction is known to lead 



to discrete collective modes inside the BCS gap [17], 
Similarly, the energy dispersion of the pairing interac- 
tion and pairing amplitude A p that becomes prominent 
at strong coupling [23| , leads to discrete collective mode 
spectrum (see below). While the exact form of this dis- 
persion is sensitive to the specifics of the strong cou- 
pling problem, it is established in the literature that, 
generally, both effects can occur near Feshbach reso- 
s S EH EL EH El- Although our understanding 
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of the detailed microscopic picture may be hampered by 
the nonpertubative nature of the strong coupling prob- 
lem, we shall see that within a simplified model used be- 
low the inequality Tito < 2 A is fulfilled under very general 
conditions. 

The dissipationless BCS dynamics 22, 25] and the pos- 
sibility to realize it in cold gases [26| attracted much at- 
tention recently [3, l2 7| . [2 8| . These investigations, with 
the exception of Ref. [24[, focused on the case of pair- 
ing interaction which is constant in the entire fermion 
energy band, concluding [29|, |3Cj that several interesting 
dynamical states, synchronized and desynchronized (or 
dephased), can be realized by a sudden change in the 
interaction strength (see the phase diagram in Ref. [29|). 

In contrast, as we shall see below, the dephased be- 
havior is suppressed in the strong coupling regime when 
due to the energy dispersion of the pairing interaction 
the Higgs mode falls inside the BCS energy gap. Un- 
der these conditions an undamped Higgs mode can be 
excited upon a sudden change in interaction. By analyz- 
ing the limit when the interaction dispersion disappears 
we show how the different regimes of Ref. [29| are recov- 
ered. This correspondence suggests an inter pret ation of 
the dephased oscillations discussed in Refs. [22l I29I l30j 
as a manifestation of the Higgs mode, algebraically de- 
phased at Twj = 2 A. 

We shall analyze the pairing dynamics in a spatially 
uniform system using the pseudospin representation [25 1 



of the BCS problem in which spin 1/2 operators s p = 
Sp ± is^ describe Cooper pairs (p, — p): 



H = - 2e " s 



p~D A pq(*)V 



(1) 



pq 



where e p is the free particle spectrum. The interaction 
Apq(t) that models the energy dispersion at strong cou- 
pling is taken in the form of a sum of a dispersing and 
nondispersing parts 



Apq(t) 
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(ai + a 2 / P /q) , / P = 
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where the dimensionless parameter g(t) specifies the in- 
teraction time-dependence, the constants > satisfy 
a\ + a,2 = 1, and vp is the density of states at the Fermi 



level. The second term in ([2]) features dispersion on the 
energy scale 7. Our motivation for choosing the model 
@ was two-fold. Firstly, the form @ is general enough 
to provide insight into the role of different features, such 
as the energy dispersion (which is controlled by the pa- 
rameter 7) and separability (which is absent unless oq 
or 0,2 vanishes). Secondly, our numerical method utilized 
the rank two form of l|2|). allowing for substantial speedup 
that could not be implemented for a more general inter- 
action A pq . In addition, the model (|2|) is physically mo- 
tivated by the theory of BCS pairing in the simultaneous 
presence of a retarded and non-retarded interaction 31 1 . 

Within the mean-field approximation, the dynamical 
equations derived from Eq.JT]) assume a Bloch form: 



dr p 

dt 



= 2b p x r p , 



bp = -(A*,A p ,e p ), 



(3) 



where r p = 2(s p ) are Bloch vectors, and the effective 
magnetic field b p depends on the pairing amplitude A p . 
The latter is defined self-consistently: 



A p = A p + zA p 



EApq(t) 
2 ' 

q 



(4) 



The interaction time dependence of interest is a step-like 
change from the initial value gi to the final value g. With- 
out loss of generality, the phase of the order parameter 
can be chosen equal zero, allowing us to consider only the 
x-component of the pairing amplitude, A p = Ap. As an 
initial state we take the paired ground state 
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The equilibrium energy-dependent amplitude A p is de- 
termined by the self-consistency equation 



Aq 



(6) 



in which A pq is given by |[5J) with the parameter values <?i 
and g for the initial and final state. The corresponding 
equilibrium pairing gap values, A p and A p , are found by 
numerically solving the integral equation (|6|). Through- 
out the paper we use the equilibrium value of the pairing 
gap at the Fermi level, A^, at the final coupling g as a 
natural energy scale to parameterize the dynamics. 

We integrate Eqs.([3]) using the Runge-Kutta method 
of the 4-th order with a time step adjusted to achieve 
sufficient precision of the calculation. In our simulation 
we use N — 10 4 , 10 5 equally spaced energy states within 
bandwidth W, —W/2 < e p < W/2, with the level spacing 
much smaller than all other energy scales in the problem. 

We analyze the quantity Ap(t) which at long times 
oscillates between the maximum and minimum values 
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FIG. 2: Long-time behavior of Af(£), the pairing amplitude 
at the Fermi level, oscillating between A+ and A_ . Shown are 
two examples of A± as a function of the initial state for a non- 
separable (circles) and a separable (squares) interaction 
Parameters used: ai — — 0.5, g — 0.43, Af/W — 0.016, 
j/W = 0.01, and oi = 0, a 2 = 1, g = 0.61, A F /W = 0.005, 
y/W = 0.01, respectively. Inset: Linear fit of a sample trace 
Af(t) vs. t" 1 / 2 used to extract A±. 

A + and A_. To find the asymptotic values A± we em- 
ploy the numerical procedure sketched in Fig. [5] inset: 
A± are obtained from the linear fits to the maxima and 
minima of Ap vs. i -1 / 2 intersection with the y-axis. 
The i -1 / 2 time parameterization is motivated by the de- 
phasing law 5A(t) oc i -1 / 2 found in Refs. 2^, 2M f° r the 
energy-independent interaction. Should the dephasing 
occur, the asymptotic values would coincide, A + = A_. 

In contrast to the above, for the interaction © the 
dephased behavior is suppressed. Instead, as illustrated 
in Fig. f2J we observe non-decaying periodic oscillations 
for a wide range of initial states, both for the initial states 
close to the normal state (<?; <C g) as well as for the 
initial states near equilibrium (t^ « g). At increasing <?j 
there is a critical point at which the asymptotic pairing 
amplitude A± becomes zero. 

To understand the origin of the oscillatory behavior for 
the dispersive interaction, Eq.©, we develop perturba- 
tion theory near the point <?.; = g. Linearizing the Bloch 
equations and taking a harmonic variation of the pairing 
amplitude, SAp ,v (t) oc e~ iart <5Ap$ , we find two collective 
modes for the x and y components of A p corresponding 
to the order parameter amplitude and phase variation 
(see Ref.[l9j]). The amplitude (Higgs) mode with fre- 
quency u) obeys the integral equation 
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(7) 



where A p is the equilibrium gap obtained from Eq.([S]). 
The equation for <5Ap" (the phase mode) is similar to 
Eq.© except for the denominator of the second frac- 
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FIG. 3: (a): The Higgs mode frequency obtained from the 
simulation with gi ~ g (circles) and from Eq.© (solid line) 
as a function of the dispersion parameter 7. The quasiparticle 
energy minimum (dashed line) lies above the collective mode 
frequency (parameters of the simulation: g = 0.61, ax = 02 = 
0.5). (b): Frequency of the Higgs mode as a function of the 
initial state for non-separable and separable interactions with 
the same parameters as in Fig[T] The frequency changes away 
from gi — g as the amplitude of oscillations increases (FigdJ , 
indicating unharmonicity of the Higgs mode. 

tion which is e 2 — cj 2 /4. As expected from Goldstone 
theorem, the equation for SA^ is solved by u> = 0. 

To find the frequency u> of the Higgs mode, we note 
that for the interaction A pq given by ©, which is an op- 
erator of rank two, Eq.([7]) turns into an algebraic equa- 
tion involving a 2 x 2 determinant. Solving it we find that 
for <22 > the frequency w lies within the BCS gap, as il- 
lustrated in FigJ3^,. To gain more insight, let us consider 
a separable interaction, a± = 0, 02 = 1, which yields 
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where A q oc / q . Balancing the factors under the sum in 
order to obtain unity on the left hand side, and noting 
that without the second factor Eq. ([5]) would be identical 
to Eq.©, it is easy to see that ui < 2Ap, i.e. the Higgs 
mode is discrete. 

Notably, as FigJS^, illustrates, the frequency obtained 
from Eq. |7J) coincides with the frequency of oscillations in 
Ap(t) obtained by simulating BCS dynamics at g « gi, 
proving that the observed excitation is indeed the Higgs 
mode. Furthermore, for g away from gi the frequency 
extracted from Ap(t) varies with g, decreasing below the 
value at g gi and approaching zero at gi g and gi 
g (see Figj3jD). This indicates unharmonicity of the Higgs 
mode that sets on at a large amplitude of oscillations. 

To test these ideas further, we considered the regime 
when the Higgs mode is strictly inside the quasiparticle 
continuum, which can be realized in the model © with 
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FIG. 4: Quenching of dephasing for weakly dispersing in- 
teraction. Asymptotic values of the pairing amplitude A± 
for a dispersing (circles) and non-dispersing (red line) inter- 
action for different initial states. The onset of dephasing is 
marked by arrows. Parameters used: a± — <J2 = 0.5, g = 0.33, 
A F /W = 0.02, j/W = 0.1, and ai = 1, a 2 = 0, g = 0.33, 
Af/W — 0.05, respectively. 



the second term of a repulsive sign, 02 < 0. In this case 
Eq.© has no real- valued solution in the region to < 2 A. 
Simulating the BCS dynamics near gi ~ g we find that 
A(t) exhibits exponentially decaying oscillations of the 
form e~ nt cos(k/t + (/>) corresponding to a complex- valued 
frequency u>. For a 2 — the collective mode frequency 
lu = 2Ap lies at the edge of the quasiparticle continuum. 
This property was linked to algebraic Landau damping 
of this mode in Refs.fH 



The discrete Higgs mode makes the BCS dynamics un- 
damped for g near gi even for weakly dispersing interac- 
tion A pq . It is interesting to connect this behavior to the 
dephased BCS dynamics found in the case of constant in- 
teraction. This is illustrated (Fig@} by the dynamics at 
weakly dispersing interaction 7 > Af , where we observe 
that the region of dephased dynamics shrinks, with the 
onset of dephasing shifting towards small g < gi. While 
the oscillation amplitude |(A + — A_) is now finite, it re- 
mains small due to dephasing in the transient region (see 
FigOJO . This behavior is consistent with the Higgs mode 
approaching the quasiparticle continuum boundary. 

In conclusion, we have shown that the energy disper- 
sion of pairing interaction leads to quenching of dephas- 
ing of the BCS dynamics, making the Higgs mode of the 
pairing amplitude discrete. Parametric control of inter- 
action in the strong coupling regime near a Feshbach res- 
onance of cold atoms can be used to excite this mode. 
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